Exponential suppression of interlayer conductivity in very anisotropic 
quasi-two-dimensional compounds in high magnetic field 
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It is shown that in rather strong magnetic field the interlayer electron conductivity is exponen- 
tially damped by the Coulomb barrier arising from the formation of polaron around each localized 
electron state. The theoretical model is developed to describe this effect, and the calculation of the 
temperature and field dependence of interlayer magnetoresistance is performed. The results obtained 
agree well with the experimental data in GaAs/AlGaAs heterostructures and in strongly anisotropic 
organic metals. The proposed theory allows to use the experiments on interlayer magnetoresistance 
to investigate the electron states, localized by magnetic field and disorder. 
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I. INTRODUCTION 

Layered materials is a wide class of compounds of very 
strong scientific and technological impact. These mate- 
rials are composed of a stack of two-dimensional (2D) 
conducting layers with weak interlayer electron coupling, 
determined by the interlayer transfer integral t^. Such 
layered structure may be of natural molecular origin, as, 
e.g., in cuprates, pnictides, organic metals or graphite, or 
of artificial origin, as in heterostructures or intercalated 
graphites. The small value of t^, which is much smaller 
than the in-plane Fermi energy Ep, provides these com- 
pounds with highly anisotropic electronic properties. 

Magnetoresistance is traditionally used to explore the 
internal electron structure of metals. [ll-[l] Magnetore- 
sistance in layered compounds has many special fea- 
tures. These are beats of magnetic quantum oscilla- 
tions (MQO),[ll angular magnetoresistance oscillations 
(AMR0),[3-@ slow oscillations of magnetoresistance 0, 
111, the phase shift of beats of MQCfS", '^l etc. The 
angular dependence of magnetoresistance is widely used 
to investigate the electronic structure of various layered 
compounds: organic metals (see, e.g., Refs. |10l - [l3j for 
reviews), cuprate high-temperature superconductors. 
\vi\ heterostructures [ll] etc. The AMRO are well 
understoodjl, @ in the framework of the standard semi- 
classical 3D theory, based on the Boltzmann transport 
equation, which gives the Shockley-Chambers formulallj 
for the angular dependence of magnetoresistance. The re- 
cent theoretical works[l3,[l^[23| on AMRO only provide 
the calculations based on this formula, which are useful 
for its practical application to extract the Fermi-surface 
parameters from the experimental data. 

There is one important feature of the angular depen- 
dence of interlayer magnetoresistance, observed in vari- 
ous strongly anisotropic layered compounds, which can- 
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not be described in the framework of the standard 3D 
theory. This is the so-called coherent-incoherent tran- 
sition of interlayer magnetoresistance, induced by dis- 
order and by magnetic field component Bz perpendicu- 
lar to the conducting lavers. |2l| - |24j This transition oc- 
curs in layered compounds with extremely high quasi- 
2D anisotropy, namely, when the interlayer transfer 
integral is less than the cyclotron energy huc — 
heBz/m*c and than the Landau level (LL) broadening 
To due to impurities. This coherent-incoherent transi- 
tion or crossover is sometimes accompanied by the metal- 
insulator transition(2ll. [25j and is always characterized 
by a strong growth of interlayer magnetoresistance R^^ 
with the increase of i3r. [2l| - [26| while the component 
of magnetic field parallel to the conducting layers only 
slightly affects the interlayer conductivity at weak mag- 
netic field. p3l[23| It is important that Rzz {Bz) grows not 
only in the maxima of MQO of magnetoresistance, when 
the chemical potential is situated in the gap between LLs, 
but also in the minima of magnetoresistance, when the 
chemical potential is exactly on the LLs and when there 
is no any gap of electron states at the Fermi level. The 
2D formulat27|| for the Shubnikov-de Haas effect, based 
on the standard 3D theory of interlayer magnetoresis- 
tance and applied to the extremely anisotropic quasi- 
2D limit, does not describe this effect. Moreover, the 
first direct comparison between the magnetoresistance in 
the coherent 3D and the so-called "weakly incoherent" 
regimes did not reveal any considerable differences. f28j 
Recently it was shown. [29 -[31 that although the initial 
model in Ref. [2^ is valid for the weakly incoherent 
limit, the calculation of the electron Green's functions 
in disordered 2D layer in Ref. [28!| is incorrect, because 
the impurity scattering is incorrectly treated via a con- 
stant imaginary part of the electron self-energy similar 
to 3D case. More rigorous calculations show (see, e.g. 
in Ref. (32j ) that in strong magnetic field, in addition 
to the oscillations due to LL quantization, the imaginary 
part of the electron self energy monotonically increases 
with magnetic field cx -v/Sr-tS^I This monotonic increase 
leads |29M3ll | to the similar increase of the Dingle temper- 
ature and of the monotonic part of magnetoresistance 
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cx \fBl^ which also changes the angular dependence of 
magnetoresistance. This improvement pol Isoj of the the- 
ory can explain some experimental observations 33j of 
the monotonic growth of interlayer magnetoresistance as 
function of B^. when V^^tz ^ huc- However, in many 
compounds the much stronger than cx \fBl growth of in- 
terlayer magnetoresistance has been observed. [21U26!] To 
explain this rapid growth of R^z (Bz), the variable-range 
hopping mechanism of interlayer electron transport has 
been proposed, [sij similar to the in-layer transport in the 
regime of quantum Hall effect (QHE). Although a rea- 
sonable agreement with experiment has been achieved 
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the underlying theoretical model is still un- 
clear. First, the actual electron interlayer hopping is not 
variable-range, as in Ref. [33 |. but always to the adja- 
cent layer, i.e. on the interlayer distance d. Moreover, the 
AMRO still exist in this limit, which means that during 
the interlayer hopping the in-plane electron wave func- 
tion is conserved, i.e. the interlayer tunnelling term in 
the Hamiltonian conserves the in-plane electron momen- 
tum. This Hamiltonian is given by Eqs. ([I])-© below. 
Without electron-electron (e-e) and electron-phonon (e- 
ph) interaction it has been studied in Refs. [29l - l3l| (see 
Eqs. (8)-(10) of Ref. 29]). However, the comparison 
with some experiments shows that the results of Refs. 
[29| - |3]| are applicable only in the intermediate magnetic 
fields, when Tq^Iz < hiUc, but the magnetic field is still 
too weak to induce the metal-insulator transition. 

The aim of the present paper is to generalize the model 
of Ref. [13 to describe the strong increase of interlayer 
magnetoresistance and the metal-insulator transition in 
very high magnetic field. 



II. THE MODEL 

The quasi-2D layered electron system with disorder 
and interaction is described by the following Hamilto- 
nian, which contains four main terms: 



Hi, 



(1) 



The first term Hq is the Hamiltonian of the noninteract- 
ing 2D electron gas in magnetic field summed over all 
layers. The second term in Eq. ((T]) gives the coherent 
electron tunnelling between two adjacent layers: 



{r) + ^]_,{r)'f,{r)], (2) 



where '^j{r) and are the creation (annihilation) 

operators of an electron on the layer j at the point r. 
This interlayer tunnelling Hamiltonian is called "coher- 
ent" because it conserves the in-layer coordinate depen- 
dence of the electron wave function (in other words, it 
conserves the in-plane electron momentum) after the in- 



Hi = J2j d^rV^ (r) ^'t(r)^'(r) 



(3) 



gives the electron interaction with impurity potential. 
The last term Hint describes the electron-electron and 
electron-phonon interaction. The presence of this term 
differs our model from that in Refs. [29l - [3l| . Even with- 
out Hint and for the case of point-like impurity potential, 
the exact solution of this Hamiltonian is not achievable. 
After inclusion of the interaction term Hint , we can only 
catch the main physical effects in this system to describe 
the interlayer electron transport. 

In the limit, <S^Tq, HuJc, the interlayer hopping tz can 
be considered as a small perturbation for the uncoupled 
stack of 2D metallic layers. The 2D disordered electron 
system in magnetic field has been extensively studied in 
connection to the quantum Hall effect. With the increase 
of magnetic field the crossover from the diffusive dynam- 
ics and classical Shubnikov-de Haas effect to the in-plane 
localization of electron states takes place. In Ref. [Ill 
this crossover is described in detail, and the value of the 
crossover field Be is estimated for the various types of im- 
purity potential. Although for point-like impurity poten- 
tial the crossover field Be was not calculated in Ref. (ssj , 
the general criteria that Be corresponds to field when 
the nearest Landau levels (LLs) become separated, i.e. 
when fkjjc ^ Tq, can be used. During this crossover the 
localization length ^ of the 2D electron states decreases 
exponentially from infinity to the value (35l - l43 | 

e « Rc, (4) 

where the classical cyclotron (Larmor) radius 

Rc ~ hkpc/eBz — kpl\j = 2v jkp^ (5) 



kp 'va the in-plane Fermi momentum, l}i = ^ hc/eBz 
is the magnetic length, and the LL filling factor v — 
F/Bz, where F = kpch/2e is the MQO frequency. In 
the quantum limit 1^ < 1 the localization length ^ ~ Ih- 
In many layered materials the in-plane electron disper- 
sion is anisotropic, and the in-plane FS has the shape of 
elongated ellipse with semi-axes given by the wave vec- 
tors kpi and kp2. Then, in the absence of impurities, the 
electron classical orbit is also an ellipse with the main 
radii 



Rci ~ 2v/kFi , Rc2 ~ 2v/kF2. 



(6) 



This equation is a particular (elliptic) case of the well- 
known property that in magnetic field the electron or- 
bits in coordinate and momentum spaces differ by the 
90 degrees rotation only. The in-plane anisotropy of the 
localization length ^ is the same as the anisotropy of the 
cyclotron radius i?c, ii cc Rd, because the electron diffu- 
sion along the conducting layer goes via the jumps from 
one cyclotron orbit to another due to impurity scatter- 
ing, and the average jump distance is proportional to the 
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size of the electron orbit in this direction. Hence, the 
locahzation length in two main directions is 



Rc 



(7) 



Any localized electron state is covered by a polaron, 
which conies from the static lattice deformation and from 
the rearrangement of the neighboring electrons. Since 
the electron localization length ^ ~ i?c is larger than the 
Debye screening radius, the charge density p (r) induced 
by this polaron follows the square of the electron wave 
function amplitude (r)|^. Hence, this polaron creates 
a positive charge distributed on the area 



(8) 



Along the interlayer direction the induced charge also 
follows the electron wave function |^e(2;)|^ and is dis- 
tributed within one crystal layer, because otherwise there 
would be a nonuniform charge distribution along the z- 
axis. This polaron lowers the energy of each localized 
electron. 




FIG. 1: The schematic view of the charge distributions (dark- 
gray areas) of polaron (lower layer) and of the electron (upper 
layer) after an electron jump to the adjacent layer. 

To jump to the adjacent conducting layer the electron 
now must overcome the Coulomb attraction of polaron, 
which creates the potential barrier. To estimate the value 
Ec of this Coulomb barrier, consider the following elec- 
trostatic problem, schematically shown in Fig. [TJ The 
flat capacitor has charge e separated by the interlayer 
distance d. The energy stored in this capacitor is 

Ec ~ 2Tre^d/eAo, (9) 

where e is the dielectric constant of the media. 



III. ESTIMATES OF INTERLAYER 
MAGNETORESISTANCE AT THE 
METAL-INSULATOR TRANSITION 

The Coulomb energy Ec, given by Eq. (0), affects the 
interlayer electron transport when it becomes comparable 
to the band- width in the interlayer direction: Ec > itz- 
Substituting Eqs. (|n])-(IHl) into Eq. (0), one can esti- 
mate the crossover magnetic field Be, above which the 



Coulomb energy hampers the interlayer electron trans- 
port: 



Be 



Fy/Set^/e^dkpik 



F2, 



(10) 



where F — kpikp2ch/2e is the MQO frequency. Some- 
times, it is convenient to use the crossover LL filling num- 
ber Vc, which corresponds to the field in Eq. ([TOl) : 



z/c ~ \/ e'^dkFikF2/8etz 



(11) 



At B > Be, or a.t < i^c, the interlayer electron trans- 
port becomes thermally activated. One could expect that 
at B > Be 



{B,)/a,, {0)^e^pi-Ee/T). 



(12) 



This simple formula does not take into account two ef- 
fects. First, the finite interlayer transfer integral effec- 
tively reduces the activation energy Ee needed for an 
electron to jump to the next layer by the bandwidth 4i^. 
This can be approximately taken into account by an ef- 
fective increase of temperature in Eq. (|12l) according to 
the rule 



T ^ T* K T + At, 



(13) 



Second, the finite probability of interlayer jumping de- 
localizes the electrons, which increases the effective in- 
plane electron localization length ^ and reduces the 
Coulomb energy Ee in Eq. ([9]). Third, the electron 
localization length vary from one state to another, and 
one actually has a distribution of localization length with 
function P (^) . This function has maximum at ^ ~ i?c 
and decrease rapidly at ^ ^ Re- However, the tails of 
the function P (^) , which play the crucial role for the 
in-plane conductivity, are also important for the inter- 
layer transport because of the exponential dependence 
on ^ in Eq. (fT^ . It can be argued that the probability 
that the electron localization area A is much larger than 
Aq « TTReiRe2 dccrcascs exponentially: 



PiA)^iA/Ao)cxpi-A/Ao). 



(14) 



Combining Eqs. and we obtain the sup- 

pression factor of interlayer conductivity 



Ozz (0) 



dA A 
— — exp 

Aq Aq 



~A 
AT 



exp 



~2'Ke^d 
eT*A 



= X^K^iX) /2 = a{X), 



(15) 



where K2 {X) is the modified Bessel's function of the 
second kind, and 



(16) 



X = ^87re2d/er*Ao = ^jAEe/T*. 
At X > 1 Eq. (US]) simplifies to 

(Tzz {Bz) lazz (0) - {X/2f'^ exp {~X) . (17) 

In Fig. [2]we plot the inverse function Rzz {Bz) / Rzz (0) — 
Ozz (0)/a,, {Bz) <x a-^ {X). 
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FIG. 2: The function q"^ (X) from Eq. GHI. The predicted 
ratio Rzz (Bz) /i?22 (0) ~ (X) is given by the same curve, 
because X (x B^, as follows from Eqs. (fT5)) .(f5)) and (0. 



where one only has to change the axis labels: on the ordi- 
nate axis {X) — >■ R^z (Sz) /Rzz (0), and on the hor- 
izontal axis X Q.iBz [T]. The interval < < IST 
of magnetic field in Fig. Id of Ref. fTsJ corresponds to 
the interval < X < 4 in Fig. [2j This calculated depen- 
dence Rzz {Bz) agrees very well with the experimental 
data plotted in Fig. Id of Ref. [18] at 6* = 0°, provided 
that there is a constant upward shift of the experimen- 
tal curve Rzz {Bz) in Ref. [l^ which, probably, came 
from in-series resistances in the experimental setup. The 
crossover field, given by Eq. ([TU|. Be « 0.5T also agrees 
with the experimental data in Fig. Id of Ref. 18|. Unfor- 
tunately, there is no experimental data on the temper- 
ature dependence of the interlayer magnetoresistance in 
Ref. [l^ to compare with out theoretical prediction. 



IV. COMPARISON WITH EXPERIMENTS 



The increase of interlayer magnetoresistance described 
above is rather general. It has been observed in many 
compounds, including heterostructures [l8| and various 
organic metals. |2l| - [26l |45| This effect, probably, also oc- 
curs in pnictides and cuprate high-Tc superconducting 
materials. For comparison with experiment we choose 
two completely different materials: the GaAs/AlGaAs 
semiconductor superlattice[18[ and the organic metal /?"- 
(BEDT-TTF)2SF5CH2CF2S03.'255 This choice is stipu- 
lated by the small value of the interlayer transfer inte- 
gral tz in both these materials, which makes possible the 
in-plane electron localization even at moderate magnetic 
field. In GaAs/AlGaAs superlattices the interlayer trans- 
fer integral tz can be controlled artificially. In the sample 
#2 of GaAs/AlGaAs in Ref. ^tz = O.OSmey = Q.ZbK 
and the interlayer distance d ^ 300i. In /3"-(BEDT- 
TTF)2SF5CH2CF2S03 the small value of tz < 0.2K ap- 
pears naturally because of the large interlayer distance 

iT.si.il 



A. GaAs/AlGaAs heterostructures 



The area electron concentration 
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11, 



= kjp/2'n- = 6 X 



in the sample #2 of GaAs/AlGaAs in Ref. [H 



corresponds to the Fermi momentum kp = 2 x 10^ cm~^. 
This gives the cyclotron radius Rc = hkpc/eBz = 
VdQnm/Bz \T]. Note that the magnetic length Ih = 
yJhcjeBz ~ 2.5 X IO^^cto/ \JBz [T] is of the same order 
as i?c, and the filling factor v ^ \. The MQO frequency 
F = k\chl2e w 11. 7r. The MQO and QHE are not ob- 
served because they are smeared completely by high tem- 
perature T = 75K. However, the localization of the elec- 
tron states is preserved and affects the electron transport. 
Substituting s = 12 A we obtain Ec ~ 27re^d/enRl « 
1.6K X B^ [T] and X = ^JAEJT* w 0.35^ [T]- Sub- 
stituting this into Eq. we obtain the dependence 
Rzz {Bz). This dependence Rzz {Bz) is shown in Fig. [21 



B. Organic metal /3-(BEDT-TTF)2SF5CH2CF2S03 



The in-plane Fermi surface in the organic compound 
;9-(BEDT-TTF)2SF5GH2CF2S03 is an elongated ellipse 
with areajl^ nkFikF2 = 0.0192 which is about 
5% of the first Brillouin zone. The interlayer distance 
d ~ 17.5^. Unfortunately, we are not aware of any 
measurements of the dielectric constant e in /3-(BEDT- 
TTF)2SF5CH2CF2S03. In organic metals the dielec- 
tric constant s varies in the interval from 8 to 50 and 
may depend on magnetic field, pressure and other ex- 
ternal parameters. [431 Therefore, for the estimates be- 
low, we take the average value e ~ 20. The exact 
value of the interlayer transfer integral tz in /3-(BEDT- 
TTF)2SF5CH2CF2S03 is not known yet, but the ex- 
perimental estimates of the maximum value of tz give 
tz < l8.5fieV = 0.21K. [13] Below we take tz = 0.2K. 
The other parameters are d = 17.5 A and F = 200 Tesla. 
Substituting this into Eqs. ^TU\\ or pT|) we obtain the 
critical magnetic field Be ~ 8.5r and the LL filling num- 
ber Uc ~ 23 in a qualitative agreement with the exper- 
imental results in Refs. [2l|, [25|. Substituting the pa- 
rameters to Eq. © we obtain Ec ~ 0.012 [K] x [T] 
and X = ^/AEJT^ « Bz [T] v/0.05 [K]/ {T + Uz). 
Substituting this into Eq. (flSl) we obtain the de- 
pendence Rzz {Bz,T) of the background magnetoresis- 
tance, averaged over MQO. This predicted dependence 
Rzz {Bz) /Rzz (0) is shown in Fig. [31 for six different val- 
ues of temperature: T = 0.66 K, 0.565 K, 0.455 K, 0.218 
K, 80 mK, and 47 mK. These values of temperature are 
taken to be the same as in Fig. 1 of Ref. [29], to sim- 
plify the comparison with experiment. As in Fig. 1 of 
Ref. [2I] , the highest curve in Fig. [3| corresponds to the 
lowest temperature. Comparison of Fig. [3land Fig. 1 of 
Ref. [15 shows, that the agreement of the above theory 
with experiment is very good, both in the temperature 
and field dependence of magnetoresistance. 
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FIG. 3: The predicted field dependence of the background 
(non-oscillating) magnetoresistance R^z (Bz) / Rzz (0) in /?- 
(BEDT-TTF)2SF5CH2CF2S03 at = 0.27^ for six values 
of temperature: T = 0.66 K, 0.565 K, 0.455 K, 0.218 K, 80 
mK, and 47 mK. The lower the temperature, the larger the 
magnetoresistance Rzz (Bz), i.e. the highest curve in this plot 
corresponds to the lowest temperature. The values of temper- 
ature are taken to be the same as in Fig. 1 of Ref. [2^ . 



V. CONCLUSION 

Above we developed a theory to describe the interlayer 
magnetoresistance in the extremely anisotropic quasi-2D 
layered metals. Our model takes into account the po- 
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